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A two-site cluster generalization of the Hubbard model in large dimensions is examined in order to 
study the role of short-range spin correlations near the metal-insulator transition (MIT) . The model 
is mapped to a two-impurity Kondo- Anderson model in a self-consistently determined bath, making 
it possible to directly address the competition between the Kondo effect and RKKY interactions in 
a lattice context. Our results indicate that the RKKY interactions lead to qualitative modifications 
of the MIT scenario even in the absence of long range antiferromagnetic ordering. 
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I. INTRODUCTION 

The competition between the Kondo effect and the 
RKKY interactions is a recurring theme in many of the 
most interesting phenomena associated with the physics 
of strong electronic correlations. When the RKKY inter- 
actions predominate, the result is long-range magnetic 
ordering, as found in many heavy- fermion materialso. In 
situations where magnetic ordering is absent, the mani- 
festations are more subtle, but often equally fundamen- 
tal. In particular, it has been suggested! 2 ] that this compe- 
tition lies at the core of the proposed "two-fluid" behav- 
ior, and "micromagnetism" found in some non-magnetic 
heavy fermion systems. Another interesting class of sys- 
tems where both the Kondo effect and the RKKY cor- 
relations are believed Jto be crucial are exemplified by 
doped semiconductors^ near the metal-ieswlator transi- 
tion. In these systems, non-Fermi liquidaEI metallic be- 
havior is observed, suggesting the coexistence of local 
moments and conduction electrons that seem decoupled 
from each other - another manifestation of the "two- 
fluid" behaviom. There are many further examples where 
these effects are of key importance. Unfortunately, there 
are very few theoretical approaches that are able to treat 
both the Kondo physics and the RKKY correlations on 
the same footing and provide a convincing picture of 
these interesting phenomena. 

Theoretically, much of our current understanding of 
strongly correlated metallic phases relies on a variety of 
mean-field descriptions, most of which essentially em- 
phasize the Kondo aspect of the problem. Several ap- 
proaches have been proposed, but the most elaborate one, 
combining many of the pre-existing ideas in the field, is 
based on taking the limit of large spatial dimensionality!! 
This method represents a generalizations of simple, but 
physically transparent mean-field ideas of Bragg and 
Williamsl, as applied to interacting electronic systems. 
In this_picturc, the electron residing on a given site is 
viewedli3 as a Kondo spin which is coupled by an ex- 
change interaction to a bath consisting of thcjcmaining 
electrons. Formally, the problem is mappecE3 onto an 
Anderson impurity model supplemented by an additional 
self-consistency condition. The Kondo resonance of the 



impurity model maps te-jbiiG heavy quasiparticle band, 
and the Mott transitiontlnlJ is driven by the vanishing 
of the Kondo temperature TKondo ~ m/m*. Besides pro- 
viding an appealing, physically transparent picture of the 
correlated state and the metal-insulator transition, the 
d = oo method provides a quantitatively accurate compu- 
tational approach valid in the entire temperature range. 
As it treats at the same level both the coherent, quasi- 
particle excitations, and the incoherent collective inelas- 
tic processes, the method even allows for a description of 
fully incoherent, non- Fermi liquid metallic stateala. 

In spite of the successes of the d = oo mean-field ap- 
proach, it remains unable to address several important 
physical questions. Since it based on a mapping on a 
single-site (impurity) model, it cannot properly account 
for the competition between the Kondo effect and the 
spin-spin correlations between neighboring sites - the ef- 
fect that we have argued is crucial in a number of phys- 
ical situations. The locality inherent in this formulation 
leads to another feature that is likely to be an artifact 
of mean-field theory: tha. "pinning" of the density of 
states at the Fermi leveO. More precisely, this effect 
can be directly traced to the momentum independence^! 
of the local self-energy, reflecting the lack of spatial corre- 
lations. In the context of strongly correlated, but weakly 
disordered systemslIZI, the pinning condition was shown 
to result in a discontinuous jump of the DC conductiv- 
ity at T = — the minimum metallic conductivity. If 
the pinning is relaxed, it is conceivable that a continuous 
behavior of the conductivity would follow, thus qualita- 
tively modifying oux picture of transport near the metal- 
insulator transition^. 

In order to address the limitations of the existing 
d = oo theory, a most straightforward approach would be 
to investigate systematic 1/d corrections resulting from 
finite dimensionality. Several different method^for per- 
forming such expansions have been proposecO'tS, but 
each of these approaches result in formidable technical 
difficulties, making it difficult to address the finite di- 
mensional effects in a simple and elegant fashion. In this 
paper, we take an alternative route: we propose to ex- 
tend the existing theories in d — oo in a way that mimics 
the most important physical effects of finite dimension- 
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ality. Given the fact that the general large-dimensions 
philosophy is based on the mapping of a lattice models 
onto appropriate impurity models, the appropriate im- 
purity model displaying the relevant physics is the two 
impurity Kondo (Anderson) modeO, which is often used 
as a simplest model for the study of the RKKY-Kondo 
competition. Using standard methodsO, we can obtain a 
lattice version of this model by self-consistently embed- 
ding it in an appropriate medium. The resulting model 
is the "minimum model" that allows us to go beyond 
the limitations imposed by the conventional d = oo ap- 
proach, without performing uncontrolled or unjustified 
approximations. 

In the rest of this paper, we define and examine this 
model, and indicate how the new features inherent to the 
RKKY-Kondo competition madifjy_.the standard d = oo 
results for the Hubbard modelliifO. Specifically, we in- 
vestigate the modifications of the Mott transition in a 
single-band Hubbard model at half filling. We conclude 
that the RKKY interactions represent a relevant pertur- 
bation, relaxing the pinning condition and qualitatively 
modifying the nature of the metal-insulator transition. 



II. THE MODEL 
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FIG. 1. Lattice structure of the doubled Bethe lattice and 
effective two impurity cluster. 

We begin our discussion by defining the model that we 
consider, and derive the corresponding self-consistency 
conditions by performing the d — » oo limit. While the 
limit of infinite dimensions does notjmpose any restric- 
tions on the lattice structure studiecO, the equations be- 
come particularl y ^w pjf and easy to derive in the case 
of a Bethe latticeOEJO. The qualitative features of the 
model will be identical as on other lattices, and the result- 
ing spectral functions are closer to the three dimensional 
situation than for example on the d = oo hypercubic 
lattice. The "minimum model" that we propose is then 
obtained by doubling the Bethe lattice (with hopping t) , 
and allowing the electrons to hop between the Bethe lat- 
tices with hopping t a t. The geometry of the resulting 



lattice is shown (for coordination number z = 3) in Fig. 
(0). 

Denoting the creation operators corresponding to the 
two Bethe lattices with a\ and b\, the Hamiltonian can 
be written as 



H = -t (4r a j<r + b L b j<?) + tab^2{at a b la + h.c), 

+ Uy^X n ai<J n ai-c! + n bia nbi-a) ~ ^^^2{n a ia + n bia ) (1) 



where U is the Coulomb potential and t is the nearest- 
neighbor hopping amplitude; t ab is the hopping ampli- 
tude between the two lattices. 

It should be stressed at this point that this model 
clearly breaks translational invariance by singling out 
pairs of sites connected by hopping elements t ab . While 
this feature appears somewhat artificial in a uniform sys- 
tem in which all neighbors are equivalent, it leads to 
a controlled and non-trivial modification of the d = oo 
limfO. In contrast to the standard single-band Hubbard 
model in infinite dimensions, in which electrons solely 
undergo temporal fluctuations, our model also allowA-fpc 
spatial fluctuations. A systematic expansion in l/oI3E2l 
includes exactly these processes and the model can there- 
fore be interpreted as including some of the effects of 
finite dimensionality. 

It is clear that the model by construction enables us 
to study nearest-neighbor spin-correlations. In physical 
terms, for t ab large, the model favors the formation of 
singlet pairs (dimers) from the "a-b" sites. Interestingly, 
this symmetry breaking is not unreasonable in disordered 
systems, where each site a has another "preferred" neigh- 
boring site 6, with which dimerization will be favored. 
This notion is at the heart of the "random singlet" or- 
dering of Bhatt and Leeci, describing the singular ther- 
modynamics of doped semiconductors. Notice, however, 
that a variety of additional interpretations is possible. In 
particular, the model may alternatively viewed as a two 
band model or as two coupled layersB3. 

As usual, the problem simplifies considerably in the 
large coordination (large dimension) limit, where a map- 
ping to an appropriate impurity model is obtained. Using 
standard methodsa, we proceed by rescaling the hopping 
amplitude t as t — > -j= (m = z — 1 is the "branching ra- 
tio" of the Bethe lattice), and taking the limit m —> oo. 
The result is an effective two-impurity Anderson model 
embedded in a self-consistently determined bath. We in- 
troduce spinors c*(t) = (a* (r), 6* (r)) and the matrix 
Green function 



G a {r-r')=( 



with 



) ~(Ta a (r)bl(r')) 
-(Tb a (T)aUr')) -(Tb a (r)bUr')) 



* T (T T c(r)ct(0)) f 



(2) 
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Notice that due to spin conservation G a = 8 aa ' G CTa v . 
The effective action can then be written in matrix form 



S ef f[c a ,c*} = -^^c; n G 1 (iw„)c o 
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where the self-consistency condition reads 



(4) 



iu> n + fi -t ab 

-t ab i^n + At 



-t 2 G(icu n ). (5) 



While solving this model for general values of the pa- 
rameters represents a highly nontrivial task, we immedi- 
ately recognize some well known limiting cases. In the 
limit t a b = 0, the two Bethe lattices decouple, and the 
model reduces tp-,tiie well known single-band Hubbard 
model in d = ootil til. At half filling, this model under- 
goes a Mott transition at U = U C2 , which is preceded, 
by a formation of a pseudo-gap and the coexistenceLL3E2l 
of a metallic and an insulating solution in the region 
U = U C1 < U < U = U^,. However, the metallic so- 



lution is lower in energyEj at T = throughout the 
coexistence region, so that U = U C2 represents a true 
zero-temperature critical point where the two solutions 
merge. 

The other easily analyzable case is the noninteract- 
ing limit U = 0. Here, a band-crossing transition takes 
place, where the density of states (DOS) at the Fermi 
level vanishes continuously and a gap opens at t a b = t. 
The origin of this transition is easy to understand: in 
the atomic limit t a b >> t the DOS reduces to the two 
(bonding and antibonding) levels at E± = ±t a b- When 
the hopping t increases, these atomic levels broaden into 
bands of width » 2t, so that the gap closes when the two 
bands overlap, at t a b = t. 



III. RESULTS 

While the limit of infinite dimensions simplifies the 
original problem considerably, solving the correspond- 
ing impurity model is still a formidable task. In the 
framework of the single-band Hubbar. d , . mod el , a num- 
ber of numerically exact techniauesliil l!3£3'Ell, as well as 
simpler approximate methodsE3H3'E3 have been used. A 
particularly -simple approach first proposed by Georges 
and KotliarEHI is based on solving the Anderson impu- 
rity model using second order—perturbation theory, fol- 
lowing Yamada and YosidaErEII. Due to the additional 
self-consistency this approach, often called the "iterated 
perturbation theory" (IPT)EJ, still has non-perturbative 
character. It is exact in both limits of U — and 
U = oo and displays a Mott metal-insulator transi- 
tion. Detailed investigations, based on other numerical 
approaches demonstratedEaEa the qualitative validity of 



most IPT predictions for the single-band Hubbard model 
in d = oo. As compared to numerically exact solutions, 
IPT requires considerably less computational effort, and 
thus represents a valuable guide to the physics of d = oo 
electrons. 
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FIG. 2. Phase diagram in the U — t a b plane. Both the 
metallic and the insulating solution are locally stable in the 
coexistence region even at T — 0, in contrast to the t a t = 
situation. The curves denoted by U c i (dashed line)and U C 2 
represent the boundaries (spinodals) of the insulating and the 
metallic solution, respectively. The two solutions merge at the 
two second order (critical) points denoted by A and B. A di- 
rect, continuous transition from metal to insulator is found 
along the critical line (to the right of point B), which is qual- 
itatively similar as for [7 = (see Figs. 3 and 4). 

In the problem that we consider in this paper, one 
has to solve a two impurity Anderson model - a task 
which is considerably more difficult than the simpler one 
impurity nie|dcl. Furthermore, numerical Monte Carlo 
approachesc3 to the two impurity Anderson model have 
proven to be largely unsuccessful at the available com- 
putational level. Alternative, exact diagonalization ap- 
proaches also appear hardly feasible in the case of the two 
impurity problem, although recent developmentsO hold 
considerable promise for the near future. Taking these 
facts into consideration, we propose to begin the investi- 
gation of the problem considered using the IPT approach 
as a useful first attempt to gain insight into the RKKY- 
Kondo competition. We note however, that in contrast to 
the earlier application of the IPT approach, in the present 
case this approximation is not exact in the U — oo limit, 
even at half filling. Still, we do not expect that these 
limitations will qualitatively modify our conclusions, es- 
pecially in view of the absence of a small energy scale 
at the first order metal-insulator transition that we find. 
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The possible instances where the limitations of the IPT 
approach could be relevant will be further discussed in 
Section IV, where we also present a quantitative estimate 
for the range of its validity. 

In the following, we will concentrate on the behavior 
at half-filling, where the Mott transition takes place at 
tab = 0, and investigate the modifications induced by 
turning on t ao ^ 0. In order to apply IPT to the present 
model, we have to compute the second order perturba- 
tion theory corrections around the non-magnetic Hartree- 
Fock solution. The second-order diagonal/off-diagonal 
self-energies in this case consist of only one diagram re- 
spectively, and are given as 



E^Cr) = -U 2 Gl y {r)Gl y {-r)Gl y {T), 



(6) 



where x,y — a,b. Since the resulting equations have 
to be solved self-consistently, the solution is obtained by 
numerical iteration until convergence is found. 
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FIG. 3. Density of states for interactions a) U — and b) 
U = 1.5 for t ab = 0.0, 0.25, 0.5, 0.75, 1.0. 

We have determined the T = phase diagram of the 
investigated model at half filling using the IPT approx- 
imation, and the results are presented in Fig. (^). At 
small values of U, as t a b is increased, the behavior is 



qualitatively the same as at U — 0. A continuous transi- 
tion takes place, at a critical value of the hopping t c ab (U) 
that is found to decrease as U increases. This behavior 
reflects the fact that the gradual band broadening due 
to the Hubbard-Mott splitting tends to close the gap. 
As an illustration we display the evolution of the DOS 
as the transition is approached in Fig. (^) at U = 
(Fig.|(a)) and U = 1.5 (Fig. |(b)). [Throughout the 
paper, all energies are measured in units of the half- 
bandwidth D = 2t] 
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a) Quasiparticle residue z — , , . — 
metallic phase as a function of the interaction U for t a b = 0.2. 
b) Density of States (-ImG(0 + )) of metallic (solid line) and 
insulating (dotted line) solutions at t a b ~ 0.2 as a function of 
the interaction U. 

For larger values of U, a coexistence region of the 
metallic and insulating solutions is found, similarly as for 
t a b = 0. The metallic solution is found for U < U C2 (t ao ); 
we note the non-monotonic dependence of U C2 (tab), which 
is first found to increase, and then to decrease as a func- 
tion of tab- Thus, as compared with t a b = 0, the addition 
of the RKKY correlation is found to extend the metal- 
lic region. More importantly, we find that the metallic 
solution disappears discontinuov sly, at U = U C2 (t a b), in 
contrast to the t a b = behaviorlijliij. The solution along 
this boundary is not characterized by a low energy scale, 
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as seen by plotting the Kondo temperature (quasiparti- 
cle weight z — 1/(1 — on one site as a function of 
U, for t ab = 0.2 (Fig. ff(a)). The density of states is 
also discontinuous at this boundary (Fig.[|(b)). We note 
the pronounced dependence of the density of states at 
the Fermi energy on the interaction U - a clear violation 
of the "pinning condition" that is observed at t a b = 0. 
This behavior is seen even more clearly by plotting the 
evolution of the metallic DOS as the boundary is ap- 
proached in Fig. (|J). Clearly, in contrast to the situa- 
tion at t a b — 0, the metallic and insulating solutions do 
not merge at U = U C2 (t ao ), so this boundary cannot be 
identified with a critical line. 

Similar behavior is obtained by examining the stability 
of the insulating solution, which is found to discontinu- 
ously disappear at U = U Cl (t a b)- As we can see from 
the phase diagram, Fig. (||), the boundaries U C1 (t a b) and 
U C2 (tab) are found to join at the critical point "B" which 
is also the end of the band crossing transition critical line. 
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FIG. 5. Density of States for t a t = 0.2 and interactions 
U= 1.0,1.8,2.6,3.4. 

We thus find that in the entire coexistence region, the 
metallic and the insulating solutions merge only at two 
points: A (t a b = 0, U = U C2 ) and B. In the rest of the 
phase diagram the two solutions are disjoint from each 
other, and the transition has a first-order character. This 
conclusion can be established even more rigorously by 
examining the local stability of each solution through- 
out the coexistence region. For this purpose, we have 
developed an approach that allows us to determine the 
stability, as described in detail in the appendix. Using 
this method, we have established that both the metallic 



and the insulating solution are locally stable, supporting 
the first-order scenario. (Note that when the procedure 
is applied in the t a b — limit, we find that in the coex- 
istence region, the insulating solution is locally unstable 
with respect to the metallic solution, in agreement with 
well established resultsBEj) . 

To obtain the location of the transition line, we have 
calculated the energies of the solutions, and determined 
the line where they coincide, as shown in Fig. (|^). Since 
the solutions merge at points A and B (see Fig. @), 
the energies of the solutions have to coincide there, and 
the (first order) transition line connects those two points. 
As an illustration, the energies of the two solutions are 
plotted for t a b — 0.2 as a function of U in Fig. (||). 
As we can see, in contrast to the t a b = findings, the 
insulating solution is lower in energy for larger values of 
U, consistent with the first-order scenario. This result 
is perhaps not surprising, as the RKKY interactions are 
generally expected to stabilize the insulating solution. 
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FIG. 6. Difference between the energies of metallic and in- 
sulating solutions for t a t = 0.2 as a function of the interaction 
U. 

In line with this first order scenario, the boundary lines 
U Cl {t a b) and U C2 (t a b) should be recognized as spinodal 
lines. An interesting question is why the two solutions 
merge at t a b = 0, i. e. why is there a critical point 
there instead of a first-order transition. The existence of 
bifurcations, i. e. critical points is usually associated with 
spontaneous breakdown of some symmetry (e.g. up-down 
symmetry in the Ising model). In the case of the single- 
band Hubbard model in d = oo the relevant symmetry 
remains yet to be discovered. 

The effects of the RKKY interactions are not limited 
to the modifications of the MIT scenario. They can also 
modify the thermodynamic behavior by affecting the dy- 
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namics of the collective spin fluctuations governing the 
finite temperature response. In order to investigate this 
aspect of the problem, we have computed the specific 
heat in the metallic and the insulating phases of our 
model. To illustrate the typical metallic behavior, we 
present results for the specific heat at U = 2.5, for three 
different values of t ab = 0, 0.2, 0.4 in Fig. (g). 

At tab = we recognize the characteristic linear spe- 
cific heat at T << TKondo ~ 0.05, corresponding to 
Fermi liquid behavior, a Schottky-like peak at T ~ 
Txondo reflecting the binding energy of the Kondo singlet, 
and insulating- like behavior at T ~ U/2 due to charge 
fluctuations (Hubbard bands). 

For t ao 7^ 0, i.e. as the RKKY interactions are intro- 
duced, a new feature appears in the intermediate energy 
range: the specific heat is enhanced at T ~ J a b, reflecting 
the emergence of additional spin fluctuations with a char- 
acteristic energy corresponding to the RKKY exchange 
interaction J ab k 4t ab /U. For t ab — 0.4 and U — 2 we 
estimate J ab ~ 0.08, giving a ("Schottky") peak in the 
specific heat at Tj ~ J/2 w 0.16, exactly where the ob- 
served enhancement occurs. 
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FIG. 7. Specific heat as a function of temperature for the 
metallic solution, at U=2.5 for t ab = 0.0, 0.2, 0.4. 

This interpretation of the observed specific heat en- 
hancement finds additional support by examining the 
corresponding behavior in the insulating regime. Typ- 
ical results are presented in Fig. (||), where the specific 
heat is plotted for U = 3.5, and t a b = 0, 0.15 and 0.3. 

In absence of RKKY interactions (t a b = 0) the specific 
heat is vanishingly small at low temperatures, reflecting 
the existence of the Mott-Hubbard gap. As in the metal- 
lic phase, the addition of RKKY interactions (t a b ^ 0) 
induces specific heat enhancements in a comparable tem- 
perature range, at Tj « 2t ab /U. We also note that the 



corresponding enhancement in the metallic state (Fig. 
0) is much more spread-out in temperature, presumably 
reflecting the scattering of these fluctuations by the cou- 
pling to the particle- hole excitations. 

IV. BEYOND PERTURBATION THEORY 

The solution of our model presented in the preced- 
ing section was based on an approximate scheme for the 
impurity problem - the-perturbation theory approach 
of Yosida and YamadaE^I (YY). While this techniques 
was utjAiaed— with impressive success in previous d = oo 
studiesMiJ'[L2l, it is important to emphasize the limita- 
tions of this approach, and identify instances where most 
important problems can be expected. 
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FIG. 8. Specific heat as a function of temperature for 
insulating solution, at U=3.5 for t a t = 0.0, 0.15, 0.3. 

When applied to single-impurity Anderson models, the 
approach of YY is generally expected to be at least qual- 
itatively correct in the entire temperature range. In this 
case, the ground state is a (local) Fermi liquid, so that 
perturbation theory converges^ and finite order correc- 
tions are sufficient. The situation is more complicated 
in two-impurity models such as the two-injmurity Kondo 
(Anderson) model. Here, a critical pointEj is found at 
half- filling, separating the RKKY and Kondo regimes. 
The emergence of this critical point has a simple phys- 
ical origin. It reflects the fact that two ground states 
with different symmetry are possible, corresponding to 
the Kondo spins being compensated either by conduc- 
tion electrons (Kondo phase) or by each other (RKKY 
phase). This criticalj-aoint, which reflects a degeneracy 
due to level crossingE2l signals a breakdown of a Fermi 
liquid description. As a result, we do not expect pertur- 
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bative approaches to be accurate in the critical region. 
Indeed, if the two-impurity problem is treated in per- 
turbation theory, the critical point is washed out into a 
smooth crossover. 

Without providing a more elaborate treatment of the 
two-impurity problem, we can at least make estimates of 
the regions in parameter space where IPT could prove in- 
sufficient. Based on the information available from stud- 
ies of the two impurity Kondo modeliia in a fixed bath, 
we expect that a critical point emerges when the RKKY 
interaction J ab is comparable to the "bare" Kondo tem- 
perature T£ ondo = T Kondo (t ab = 0) 
Kondo temperature vanishes 



Since near U C2 the 
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Kondo 



in,, 
m* 



(U C2 - U), 

but the exchange interaction remains finite 



Jab ~ t 2 ab /U 



(J) 



(8) 



one can expect that increasing U at finite t ao drives the 
system from a Kondo to an RKKY metallic phase. The 
critical line where this could take place can be estimated 
by equating T^ ondo and J a b, and using Eqs. (0)-(§), we 
find 



U RK KY(tab) ~ U C2 - At 2 ab /U c 



(9) 



This expression is valid only in the t a b — » limit, 
where to leading order we have ignored the modifications 
of the (self-consistently adapting) electronic bath. This 
estimate is plotted in Fig. (|^), where it is compared 
with the perturbation theory predictions for the metal- 
lic phase boundary, and the location of the first-order 
metal-insulator transition. 
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FIG. 9. Phase diagram indicating the critical line corre- 
sponding to the Kondo-RKKY phase transition. Also shown 
is the first order line found in our model. 



As we can see, according to IPT, the first order transi- 
tion 'preempts the approach to the RKKY-Kondo critical 
line, supporting the validity of IPT-based predictions. 

The ITP prediction that the introduction of RKKY 
interactions induces a first-order metal insulator tran- 
sition finds additional support if we recall that a sim- 
ilar conclusion was obtained by introducing additional 
RKKY interactions in the large-N approaches to corre- 
lated electronsEil. However, we emphasize that this ap- 
proach did not have the two impurity Kondo physics 
build in, and, in particular, the possibility that the 
RKKY-Kondo competition induces a nontrivial critical 
point even on the impurity level. 

An interesting question that deserves further study is 
the role of the two impurity Kondo model critical point 
in the destruction of the metallic phase. Of course, this 
question would be particularly relevant if additions of 
small perturbations, perhaps disorder, could stabilize the 
metallic phase to larger values of U, so that the relevant 
critical point becomes physically accessible. 



V. CONCLUSIONS 

In this paper, we addressed the role of short-ranged 
magnetic correlations in determining the behavior of 
strongly correlated electronic systems. To account for 
these effect, which are not properly treated by existing 
approaches, we propose a two-site cluster generalization 
of the Hubbard model in infinite dimensions as the sim- 
plest model containing the relevant physics. The model 
is mapped onto a two-impurity Kondo- Anderson model 
in a self-consistently determined bath, making it possible 
to directly address the competition between the Kondo 
effect and RKKY interactions in a lattice context. 

Using a well known approximation scheme for solving 
the self-consistency conditions, we have determined the 
phase diagram of our model and discussed the modifica- 
tions of the metallic behavior. We find that the addi- 
tion of RKKY interactions induces a first-order metal- 
insulator transition, by energetically favoring the insu- 
lating phase. Additional low-energy spin fluctuations 
emerge, leading to enhancements of the specific heat in 
the intermediate temperature range, both in the metallic 
and insulating phases. 



APPENDIX: STABILITY OF d 
SOLUTIONS 



oo 



In this appendix, we describe a method that can be 
used to examine the local stability of the d = oo self- 
consistency equations, Eqs. @-(H). The method is based 
on the observation that these equations can be derived 
using a variational approach, i.e. by extremizing a certain 
functional in analogy with standard Landau-Ginsburg 
formulations of mean-field theory 

In the present case, this functional takes a form 
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(^„)], 



(10) 



where i 7 i mp [G (T (iw„)] is the free energy of the two- 
impurity model as defined by the action of Eqs. (Q) and 



F imp [G a {iu n )] = -~ln y ^ c *C ce - s ="[ G ^^")]. (11) 

Here, we consider Eqs. (f§)-(|]) as a definition of the ef- 
fective action, so that F[G a (iui n )] is a functional of a 
arbitrary, yet unspecified function Go-(iu> n ). 

In analogy with conventional Landau- Ginsburg formu- 
lation, we think of G a (iLu n ) as an order parameter (func- 
tion). The mean- field equations are then obtained by 
cxtrcmizing the above functional with respect to varia- 
tions in the form of G cr (zw rl ). The extremum condition 
reads 



SF[G] 

5G a (iu n ] 



= 0, 



giving 



G a (iu n ) = G™*(iu n ). 



(12) 



(13) 



Here, G l ™ p {iL0 n ) is the local Green function of the impu- 
rity model corresponding to a fixed "bath" Green func- 
tion G a (iu) n ) 

G™V(lUJ n ) = (4(t Wn ) c(iWn)>S.„[GI- (14) 

Note that G™ p (iw„)is also a functional of G a (iLu n ). 

Obviously, Eq. ( p^ ) is identical to Eq. (||), so that 
we recover the d — oo the self-consistency conditions Eq. 

Before going further, it is worthwhile to comment on 
the physical interpretation of the functional of Eqs. (|Io|)- 
(pd|). Using the expressions for the free energy of d = oo 
modelsa, one can show that our functional reduces to 
the free energy, when evaluated for the value of G a {iw n ) 
corresponding to the solutions of the self-consistency con- 
ditions, Eqs. (Q)-(|5|). We conclude that F [G] represents 
a free energy functional, in the usual Landau-Ginsburg 
sense. 

When the self-consistency conditions are solved numer- 
ically, one typically makes an initial guess for G a (iu n ), 
defining the effective action of the impurity model us- 
ing Eq. (^|). The impurity model is then solved by any 
method available, and a new value of G CT (ia; n ) obtained 
from Eq. (^). In numerical analysis, this iterative pro- 
cedure is known as the "substitution-iteration method" , 
which is repeated until convergence is reached. In the fol- 
lowing, we prove a general theorem that such an iterative 
procedure converges towards a local (nearest) minimum 
of the free energy functional. We note that the set of 
all possible functions G a (iuj n ) form a vector space (more 



precisely an infinite dimensional Hilbert space), and for 
notational simplicity, denote these vectors by x. 
We define the gradient vector field g(x) as 



g (x) = a x F[x], 

so that Eq. ( |l2| ) can be written as 

g(x)| x =x =0. 



(15) 



(16) 



Here, x corresponds to the solution of the self- 
consistency condition, i.e. is a locally stationary point 
of F[x]. If wc define further the quantity 



f = x g(x), 
we find that at the stationary point 
x = f(x ). 



(17) 



(18) 



In this language, the "substitution-iteration' search for 
the solution can be written as 



x(n + 1) = f(x(n)), 
and the solution corresponds to 

x D = lim x(n). 



(19) 



(20) 



Note that the increment of x(n) can be also written as 

Ax(ra) = x(n + 1) - x(n) = -g(x(n)). (21) 

As we can see, the iteration takes the vector x in the 
direction opposite to the gradient, i. e. "down the hill" 
, so that the iteration converges only in the vicinity of a 
local minimum. 

On general grounds, we expect the physical solutions 
near first-order transitions to be locally stable. We can 
check this stability, by making a small modification in 
the initial conditions that produce the respective solu- 
tions. More precisely, we should first find the convergent 
metallic and insulating solutions and x^. We can 
then examine the stability of, for example, the insulating 
solution by re-starting the iteration search from a new 
initial guess 



x l (n = 0) = (l-c)x^ + cx 



A I 



(22) 



The solution is locally stable if for c sufficiently small 
the iteration procedure converges to x = x^. We can 
similarly check the stability of the metallic solution by 
choosing c»l, 

In order to apply these ideas, we first test them in the 
well examined limit corresponding to t a b = 0. We find 
that the metallic solution is stable throughout the coex- 
istence region, but that the insulating solution becomes 
unstable as T — > 0. These findings are in complete agree- 
ment with the resulted that C/ c2 is a T = critical point 
at which (upon reducing U) the insulating solution be- 
comes unstable, and a new metallic solution emerges. We 
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can apply these ideas for t a b Oj in which case both so- 
lutions are found to be locally stable, in agreement with 
a first-order scenario. 




0.4 0.6 
Admixture of metallic solution 

FIG. 10. Landau free energy for tab — 0, U = 2.6 and 
temperatures T = 0.003 - 0.033 in steps of AT = 0.006. The 
inset shows the free energy difference between the two minima 
as a function of temperature. 

We conclude this discussion by an explicit construction 
of the "Landau-Ginsburg" functional, which represents a 
nice illustration of the above stability considerations. To 
do this, we note that in the case of the d = oo equations, 
the gradient vector takes the form 



g = G imp [G] G, 



(23) 



which can be calculated by any method that solves the 
Anderson impurity model, e.g. the YY approached. Once 
the gradient is available, it is possible to determine the 
evolution of F[x] along any particular direction in the x 
space. In particular, we expect the physical solutions to 
be local minima, separated by an unstable solution (local 
maximum or a saddle-point). It is thus useful to consider 
the direction (vector) connecting the two solutions, which 
can be parameterized as 



x(f) = (l-^+^ 



M 



(24) 



The increment of f[x(-£)] can be expressed as a line in- 
tegral 



AF(£) = F[x(£)] - F[xJ] = f 

Jo 



d£-g(x(t)). (25) 



We can numerically compute this line integral by an ap- 
propriate discretization procedure, and typical results in 
the t a b = limit are plotted in Fig. jl0|). 



Here, we show AF(£) for U — 2.6 and for several dif- 
ferent temperatures. As we can see, at T / both the 
metallic and the insulating solution are locally stable, but 
the insulating one becomes unstable as T —> 0, in agree- 
ment with our stability considerations. We also note that 
as the temperature is increased, the free energy of the 
metallic solution (I = 1) increases, until the spinodal is 
reached, where the local minimum becomes an inflection 
point and becomes even locally unstable. Of course, this 
instability is preempted by a first order transition, which 
in this case happens at , . fin ite temperature, in agreement 
with findings of RefsMlii 
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